Expansion of a wave-packet in lattices with disorder and nonlinearity 
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We show, theoretically and experimentally, the counterintuitive result that an increase of disorder 
can result in an enhanced spreading of an initially localized excitation. Moreover, we find that adding 
a focusing nonlinearity facilitates the expansion of the wave-packet even further by increasing its 
effective size. We find a clear transition between between the regions of enhanced spreading (weak 
disorder) and localization (strong localization) described by a "diffusion peak." 
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The simultaneous presence of disorder and nonlinear- 
ity in wave systems can give rise to a variety of rich and 
complex behavior [![. For a disordered system, it is a 
widespread belief that an increase in disorder will fa- 
cilitate localization, that is, the initial spreading of an 
evolving wave-packet eventually stops (HH)]. When - in 
addition to disorder - a weak nonlinearity is added to the 
system, it has been predicted that the nonlinearity will 
cither promote or inhibit the wave-packet spreading, de- 
pending on the systems details and the relative strength 
of disorder and nonlinearity Thus, a general picture 
of the effects of the interplay between disorder and non- 
linearity, remains unclear. Disorder-induced localization 
is based on wave interference and, hence, it is a univer- 
sal concept applicable to a variety of physical systems [B| , 
such as the transport of acoustic waves, microwaves, spin 
waves and matter waves In this respect, opti- 

cal waveguide lattices [l(J have emerged as ideal sys- 
tems in which the interplay of disorder and nonlinearity 
can be observed by means of simple table-top experi- 
ments. After a first experimental study of wave evolu- 
tion in a weakly-disordered nonlinear 2D fiber array [ll[ , 
a ground-breaking work 12j was published where the av- 
eraging over multiple individual realizations of disorder 
played a key role, resulting for the first time in the exper- 
imental demonstration of genuine Anderson-localization. 
In that work it was also shown that focusing nonlinear- 
ity indeed facilitates localization. Whereas this work was 
performed in a two-dimensional (2D) system, in a subse- 
quent experiment [l3| the impact of nonlinearity on local- 
ization in a one-dimensional (ID) system was experimen- 
tally analyzed, culminating in the same result: focusing 
nonlinearity indeed enhances localization of propagating 
waves. 

In this work we present theoretical and experimental 
evidence that the presence of weak disorder can promote 
the spreading of a wave-packet. Moreover, this effect is 
further enhanced by the addition of focusing nonlinearity. 
We focus on finite lattices, where the system dimensions 
can be smaller than the localization length, and monitor 



the spreading of the wave-packet during the first stages 
of its evolution. 

Longitudinal light propagation in weakly-coupled, non- 
linear optical waveguide arrays can be modeled by a set of 
normalized, discrete nonlinear Schrodinger (DNLS) equa- 
tions [HEd]: 



.du n m-^ 



(1) 



Here, u n is the amplitude of the waveguide mode in the 
n-th waveguide. The coordinate n depends on the di- 
mension and lattice type [for example, n = (k, I) in a 
2D square lattice]. The quantity e n is the propagation 
constant (i.e., "site energy" in the quantum mechanical 
context) of the n-th guide. The hopping between adja- 
cent lattice sites n and m is described by the coupling 
constant V a , m . In our model we impose disorder on both, 
the propagation constant [e n G {— W c /2, W e /2}] and the 
coupling constants [Vn,m £ 1 + {—W c /2, W c /2}}; the dis- 
order is therefore characterized by the disorder strengths 
W e and W c , respectively. Note that we limit W c to the 
interval {0, 2}, to insure that the inter-site coupling is al- 
ways positive. One of the most useful measures used for 
the description of localization phenomena in finite sys- 
tems is the participation ratio (PR) [IH which gives a 
rough estimate of the number of sites that are occupied 
by the wave-packet, i.e., where the light has a signifi- 
cant amplitude at the output facet of the sample. As 
pointed out in [la ], in order to get meaningful data for 
finite lattices, one has to also average over different re- 
alizations of disorder, and thus, work with an averaged 
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PR: R = (£n=i |u„| 4 /P 2 ^> [Hi, wherc p is thc total 
optical power, and N is the total number of waveguides. 

We start our analysis by carrying out extended nu- 
merical simulations in ID lattices. We integrate Eq. ([1}, 
for many random realizations, and perform an averaging 
process for the relevant quantities. In the ID case, we 
use N = 81 waveguides with an initially localized exci- 
tation: u n (0) = \/P<5n,n c , where n c corresponds to the 
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Figure 1. (color Online), (a) Sketch of the light propagation 
in a ID waveguide array, when only one waveguide is excited. 
(b)-(c) Simulated intensity distribution in the lattice sites for 
ordered (b) and on-site disorder level W t = 0.72 (d) after 
propagation through the waveguide array when only the cen- 
tral waveguide is excited, (d) Simulated averaged PR R vs. 
on-site disorder for different levels of inter-site disorder after 
linear (low power) propagation, (e) Simulated averaged PR 
vs. Power P and onsite-disorder W ( after nonlinear propaga- 
tion (thick contour denotes R = Ro). In all cases, the PR is 
normalized to Ro = R(W E = 0). 



input position. The participation ratio was taken at the 
propagation distance z max = 20 (in normalized units). 
In order to simulate linear propagation, the total opti- 
cal power P is chosen very small (< 0.01). Figure [Tf a) 
shows light propagation in an ordered waveguide array. 
In order to reproduce Anderson's results Q, we show in 
Fig. Qlb)-(c) output distributions for the ordered lattice 
and one level of pure on-site disorder, showing a ballis- 
tic spreading for the ordered (W e = 0) case [Fig. H^b)]. 
The diffraction pattern exhibits distinct side lobes and 
small amplitudes around the initially excited central site. 
When weak disorder is introduced [W e = 0.72, Fig.QJc)], 
the power contained in these lobes is redistributed to the 
waveguides close to the center. By increasing the disorder 
further, a localization around the excited site is observed, 



with exponentially decreasing amplitudes on both sides. 
However, an important aspect here is the occupation of 
the individual lattice sites, that defines the PR, i.e. the 
effective width of the diffraction pattern. Our simula- 
tions reveal something very surprising: As the (weak) 
disorder strength is increased, the average PR increases 
[black solid line in Fig. [ljd)]. Upon further increment 
of the disorder, the spatial profile reduces its expansion 
and the PR decreases. The figure shows a clear maximum 
expansion which separates regimes of expansion and lo- 
calization; i.e, a "diffusion-peak". We note in passing 
that, even though this phenomenon was encountered be- 
fore [IH, it was neither discussed in depth nor explained. 

We also investigate the role of mixed disorder by in- 
creasing the amount of inter-site disorder W c . The results 
for W c = W e /3, We/2, and W £ , are shown in Fig. QJd) as 
the red dashed line, the blue dotted line, and the pink 
dash-dotted line, respectively. From our simulations it is 
evident that, as the inter-site disorder increases, the ini- 
tial growth of the PR is reduced, and eventually vanishes 
(i.e. R/Rq < 1 for all disorder levels) for sufficiently high 
inter-site disorder. 

Next, we analyze the impact of nonlinearity, i.e. the 
propagation of high-power optical beams. In general, one 
would expect that a focusing nonlinear term in Eq. ([1]) 
would facilitate the self-focusing of the excitation around 
the initially excited site fl9j | . We find that, in contrast, 
in the weak-disorder regime the increase of the PR is 
enhanced in the presence of a small amount of nonlin- 
earity. Figure DJc) shows the average participation ratio 
of the profile after propagating a certain distance in the 
presence of on-site disorder W c and power (nonlinearity) 
P in a ID lattice. Interestingly, when the disorder is 
switched off and only nonlinearity remains, the picture 
is similar: the PR increases up to some maximum value 
due to redistribution of the power in the waveguides and 
then drops as the wave-packet localizes [2(|. Therefore, 
below a critical value (P < 4), nonlinearity can facilitate 
the delocalization process. 

These theoretical predictions are confirmed by our ex- 
periments, whose results are summarized in Fig. [2] The 
experiments are carried out in waveguide a rray s fabri- 
cated by the direct-writing laser technology |23| in pol- 
ished bulk fused silica wafers. The dimensions of each 
guide are 4x12 /iin 2 with a refractive index increase of 
«5x 10 -4 and a propagation length of 100 mm. For the 
analysis of ID samples, we fabricated several arrays with 
N = 81 sites each: One without disorder, and nine with 
varying degrees of disorder (and 30 realizations for each 
degree of disorder) . Disorder was created by varying the 
spacing between the guide centers: d = (23 ± r5d)[/im], 
5 d = (0,0.25,0.5,0.75,1,1.5,2,3,4,6). The difference 
in overlap of the individual waveguide modes addition- 
ally creates a statistic detuning of the guides and, there- 
fore, creates an additional on-site disorder. At the input 
facet, light at A = 800 nm from a ThSapphire laser was 
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Figure 2. (color Online), (a) Microscope images of an or- 
dered and disordered ID array, (b) Experimental intensity 
output patterns. Upper row: ordered array; middle row: sin- 
gle realization disordered array; bottom row: averaged disor- 
dered output pattern, (c) Experimental PR R versus disorder 
strengths for different levels of input power. 



launched into a single waveguide in the center of the ar- 
ray using a standard microscope objective. At the end 
facet of the sample, the intensity patterns were recorded 
with a CCD camera. Whereas the linear regime (Pu n ) 
was measured by using the continuous wave-mode of the 
laser (i.e., the laser power was only a few mW), the non- 
linear regime was studied in the pulsed regime with pulse 
peak power of P rdl ~ 1MW and P nl2 ~ 2.7MW. The 
participation number R and its standard deviation were 
computed for all the different settings, and normalized 
to the value of the linear single-site excitation in an or- 
dered array (-Rrj,Zm)- A microscope image of the front 
facet of an ordered and a disordered ID lattice is shown 
in Fig. [H[a). Fig. |^b) shows different propagation pat- 
terns of an ordered array (upper row), a single disordered 
realization (middle row) , and an averaged output profile 
for the same disorder level (lower row). From each real- 
ization, the PR was computed and then averaged. As it 
can be clearly seen [see Fig. [U[c)] , in the linear regime, 
we observe the derealization tendency with the increase 
of weak disorder and a diffusion-peak at a disorder level 
of Sd = 1.5 yum spacing variation between the guides. 
When nonlinearity comes into play, the diffusion-peak 
shifts towards smaller disorder levels, and its height in- 
creases significantly. However, the derealization is only 



observed for sufficiently small nonlinearities. If the non- 
linearity is too high, the derealization effect vanishes as 
the self-focusing is sufficiently strong to inhibit any dif- 
fusion/expansion process. 

In 2D waveguide arrays, one of the simplest struc- 
tures is the square lattice where each site is coupled to 
four nearest-neighbors only. We fabricated lattices with 
21x21 sites with three different levels of disorder and one 
ordered array. The mean distances dh or = (17 ± 8h)[^m] 
with 8h = (0,2,4,6). Although the vertical distance was 
kept fixed at d v = 23 \p,m], this setting is equivalent to 



a 2D fully disordered lattice [24|. Figure EJa) shows a 
microscope image of an ordered square waveguide lattice, 
while Fig. [3](b) shows its experimental output intensity 
pattern for a single-guide excitation. It exhibits four dis- 
tinct side lobes. Figure |3|c) shows numerical simulation 
results for the disordered lattice case. They are qualita- 
tively similar to the ones obtained in ID lattices, with 
the height and position of the diffusion peak decreasing 
as mixed disorder increases. Our experimental results, 
summarized in Fig. EJd) , validate this scenario: For small 
inter-site disorder the PR significantly grows, and for fur- 
ther increase of the disorder the PR drops, resulting in 
localization of the wave-packet (black solid line). Im- 
portantly, distinct side lobes in the diffraction pattern 
only occur when a single waveguide is excited [H, fToj ) . 
Numerical simulations show that a wide Gaussian beam 
creates a propagation pattern with a single lobe at the 
center, whose height (width) decreases (increases) con- 
tinuously during propagation. Its PR decreases mono- 
tonically with increasing disorder. This was experimen- 
tally verified by using a Gaussian input beam that cov- 
ered approximately 9 sites. We clearly see that there is 
no derealization enhancement for such initial condition 
[Fig. [3|d), red dashed line]. As for pure inter-site dis- 
order, no diffusion-peak is observed in our simulations, 
we therefore conclude that the variation of the inter-site 
spacing results in a mixed disorder in the fabricated sam- 
ples. 

Most explanations concerning diffusion processes in 
disordered lattices resort to band structure changes. 
Such is the case of disordered quasicrystal structures [16| ■ 
However, in this case the reduction of the pseudo-gaps 
due to disorder does not necessarily imply that a par- 
ticular excitation located at any input position will be 
able to excite more states. For weak disorder it is not 
possible to claim that, by considering a single-site exci- 
tation, a particular state will get excited; this may be 
only possible for strong disorder when all states are very 
localized. Therefore, in our view it is not useful to discuss 
a precise excitation of states, as this would be certainly 
a matter of probability [l5| and it will not correspond 
to a representative case. Instead, we rely on simple dy- 
namical arguments. A wave-packet initially localized in 
a single waveguide can be decomposed into a superposi- 
tion of plane waves, each one with a different transverse 
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Figure 3. (color Online), (a) Microscope image of the fab- 
ricated ordered square lattice, (b) Corresponding averaged 
output intensity profile, (c) Numerical simulation of the PR 
for different mixed disorder widths: W c = (black solid line) , 
W c = W e / 4 (red dashed line), W c = W c /2 (blue dotted line), 
(d) Experimental results for single-site (black solid line) and 
gaussian (red dashed line) input excitations, averaged over 25 
realizations. All data is normalized to the ordered array Rq. 



velocity. In the absence of disorder, during propagation 
this wave-packet will form the typical side lobes of the 
diffraction pattern 0, [uj where the waves possess a high 
transverse velocity, and most of the propagating power 
is concentrated in these lobes. For weak disorder, one 
can assume that launching light into a single waveguide 
excites a number of different linear modes with differ- 
ent transverse velocities. The main lobes still propagate 
without strong distortions as those plane waves possess 
sufficient kinetic energy to overcome the disordered po- 
tential wells and to move across the lattice. However, 
some waves with less kinetic energy will get trapped in 
the impurity regions of the particular disorder distribu- 
tion. As a consequence, the wave will expand more ho- 
mogeneously, localizing energy in the input region (small 
velocity waves), in the lobes (high velocity waves), and 
in between (intermediate velocities). When the disorder 
is further increased, the waves in the side lobes will also 
get trapped as their energy becomes comparable to the 
deeper disordered potential wells. Therefore, light re- 
duces its expansion and the localization starts to show 
its signature [see Fig. Q] and Refs. 11 -HI]. In this light, 
the PR will increase at small disorder strength and de- 
crease only at sufficiently high disorder. It is important 
to clarify that this process does not have implications in 
the edge-to-edge-diameter of the wave-packet; it rather 
implies a denser distribution of the light and therefore 
a higher PR. The side lobes remain the dominating fea- 
ture of the propagation pattern even in the far field [5l| . 
Also, the effect is not connected to the initial diffusive 
spreading at short propagation distances p2| . 



In the absence of disorder, the side lobe-dominated en- 
ergy distribution mechanism is also found in higher di- 
mensional lattices: 2D square and honeycomb, and 3D 
lattices - all of these structures will therefore exhibit the 
initial increase of the PR for small disorder and, there- 
fore, a diffusion-peak. An exception is the triangular 
(hexagonal) lattice which, due to its high coordination 
number (6), possesses an unusual discrete diffraction pat- 
tern with no distinct lobes, resembling the diffraction 
pattern for a continuous medium. Thus, our numerical 
simulation of the triangular lattice finds no enhancement 
of the PR, even at weak disorder, in agreement with the 
experimental findings in [12 [ . 

For a disorder level in the vicinity of the diffusion- 
peak, the addition of a weak nonlinearity will increase 
the (random) refractive index at each site, in an amount 
proportional to the light intensity on the site. This causes 
the high- intensity side lobes to be affected the most, 
whereas the 'slow' waves in the center are not affected 
much. However, the latter are already localized, as their 
energy content is insufficient to overcome the disordered 
potential wells. Thus, the nonlinear deepening of the ran- 
dom potential wells, renormalizes the disorder strength 
of the individual lattice sites. This, in turn, causes an 
increment of the participation ratio which occurs now 
at smaller values of disorder strength than in the linear 
case, in agreement with our experimental and numerical 
results. 

In summary, we have shown both, theoretically and ex- 
perimentally that the presence of weak disorder can lead 
to an enhancement of the spreading of an initially local- 
ized optical beam in ID and 2D waveguide arrays. More- 
over, the addition of focusing nonlinearity facilitates the 
expansion even further. The regions separating enhanced 
expansion (weak disorder) and localization (strong dis- 
order) can be clearly separated by a "diffusion peak", 
whose general behavior can be explained by dynamical 
arguments. 
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